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ABSTRACT 

The presence on the Virgo cluster outskirts of spiral galaxies with gas deficiencies as strong as those of the inner 
galaxies stripped by the intracluster medium has led us to explore the possibility that some of these peripheral 
objects are not newcomers. A dynamical model for the collapse and rebound of spherical shells under the point 
mass and radial flow approximations has been developed to account for the amplitude of the motions in the Virgo 
I cluster (VIC) region. According to our analysis, it is not unfeasible that galaxies far from the cluster, including 
those in a gas-deficient group well to its background, went through its core a few Gyr ago. The implications would 
be: (1) that the majority of the Hi-deficient spirals in the VIC region might have been deprived of their neutral 
hydrogen by interactions with the hot intracluster medium; and (2) that objects spending a long time outside the 
cluster cores might keep the gas deficient status without altering their morphology. 
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1. AN INTRIGUING POSSIBILITY WORTH EXPLORING 

A recent characterization of the large-scale 3D distribution of 
the neutral gas (H I) deficiency around the Virgo I Cluster (VIC) 
region by Solanes et al. (2002, hereafter Paper I) has shown 
that there are a significant number of galaxies with a dearth 
of atomic hydrogen at large Virgocentric distances. These pe- 
ripheral gas-deficient objects, which can be observed both in 
the cluster front and in a probable background group well be- 
hind the cluster core, show gaseous deficiencies comparable in 
strength to those measured in the centers of Virgo and other rich 
galaxy clusters. 

One of the mechanisms that can most naturally account for 
the observed reduction in the interstellar gas content of cluster 
galaxies is the ram pressure ablation caused by the rapid mo- 
tion of galaxies through the dense intracluster medium (ICM). 
There is now compelling evidence for the decisive participation 
of this process in the gaseous deficiencies of spirals located in 
the centers of rich clusters, either directly from observations 
(Giovanelli & Haynes 1985; Gavazzi & Jaffe 1987; Dickey & 
Gavazzi 1991) — including the discovery of shrunken gaseous 
disks (Cayatte et al. 1994; Bravo-Alfaro et al. 2000) and the 
finding that Hi deficient spirals are on very eccentric orbits 
(Solanes et al. 2001) — or from theoretical studies that have 
checked the efficiency of this mechanism (Stevens, Acreman, 
& Ponman 1999; Quilis, Moore, & Bower 2000; Vollmer et al. 
2001, to name only a few). 

In the outer cluster regions the low density of the intergalactic 
medium calls, in principle, for alternative gas removal mecha- 
nisms, such as gravitational tidal interactions. Indeed, observa- 
tional evidence suggests that processes of this kind might have 
played an important role on the evolution of the galactic popula- 
tion in distant clusters (e.g., van Dokkum et al. 1999) due to fa- 
vorable conditions for frequent low-relative velocity encounters 



among the galaxies in early epochs. Although in the VIC re- 
gion low-relative-velocity galaxy-galaxy interactions may also 
be responsible for the gaseous deficiencies observed in some 
of the peripheral galaxies, it should not be forgotten that the 
dynamics of the Virgo region is dominated by large-scale non- 
Hubble radial streaming motions. In this context, it is plausible 
that some galaxies at large Virgocentric distances are on very 
eccentric orbits that carry them right through the cluster cen- 
ter with high relative velocities and are therefore liable to have 
suffered a strong interaction with the ICM. 

One of the most influential studies of the Local Supercluster 
based on dynamical model calculations is the analysis by Tully 
& Shay a (1984) of the infall of galaxies in the Virgo Southern 
Extension (or Virgo II cloud) toward the VIC. The lumpy dis- 
tribution of galaxies in space led these authors to predict a very 
irregular infall rate which would be responsible for the secular 
evolution of the mix of morphological types in the cluster. It 
was suggested that the formation of the cluster took place at an 
early epoch — when the universe was about one fourth of its 
present age (R. B. Tully 2002, private communication) — by a 
first generation of, probably, early-type galaxies. Afterwards, 
infall was reduced until very recently when the large spiral-rich 
Virgo II cloud has begun to fall into the cluster diluting the 
fraction of early-type systems. The fact that Tully & Shaya saw 
very few outwardly moving galaxies outside the 6° VIC circle 
supported their argument that most, or perhaps all, spirals and 
irregulars in Virgo, mostly supplied by the Virgo II cloud, were 
recent arrivals. 

All the findings of Tully & Shaya (1984) were based on a data 
set that contained a limited number of galaxies with relatively 
uncertain distance estimates. Now, with a much larger sample 
and more accurate distances, we provide evidence that a sub- 
stantial number of galaxies, with a wide range of clustercentric 
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distances, are expanding away from Virgo. This suggests that 
these galaxies are probably reemerging after infall. If our 




FIG. 1. — Voxel projection of the 3D distribution of HI deficiency in the 
VIC region. The plot is in rectangular equatorial coordinates with distances 
given in Mpc. The xy-plane corresponds to Decl. = 0°, the x- and y-axis point 
to R.A. = 12 and 18 hr, respectively, and the z-axis points to the north. The 
central dark spot is associated with the cluster, with M87 being right at its 
center. The other two enhancements are peripheral regions of H I deficiency 
in the frontside and backside of the VIC. Our position is at the origin of the 
coordinate system. 

interpretation of the situation is correct, then there is a more 
continuous influx of galaxies into Virgo than previously antic- 
ipated, so the hypothesis that the strong gaseous deficiencies 
currently seen in some peripheral objects were originated in a 
previous infall episode is worth exploring. 

2. THE MAIN H I DEFICIENCY ENHANCEMENTS OF THE 
VIRGO CLUSTER 

A continuous representation in rectangular equatorial coor- 
dinates of the spatial distribution of Hi deficiency in the VIC 
region is shown in Figure 1 . This image is like a radiograph in 
which the shade intensity informs on the average Hi deficiency 
of the galaxy distribution observed from a given viewing angle. 
Three dark spots indicating the accumulation of galaxies with 
a dearth of neutral hydrogen are easily identified aligned along 
the line-of-sight (LOS). In Paper I it is demonstrated that the 
central enhancement of the Hi deficiency, which is essentially 
coincident with the cluster core, arises from numerous gas-poor 
galaxies with LOS distances ranging from ~ 16 to 22 Mpc. An- 
other region with important Hi deficiency is associated with 
several nearby galaxies at LOS distances <~ 10-15 Mpc, prefer- 
entially located to the north of M87, and moving away from the 
cluster with large relative velocities. Finally, there is a tentative 
background group of galaxies at LOS distances ~ 25-30 Mpc, 
most with systemic velocities close to the cluster mean, lying in 
the region dominated by the southern edge of the classical M49 
subcluster, and clouds W and W. 

While the frontside enhancement of the Hi deficiency is pro- 
duced by gas-poor galaxies that appear relatively clustered in 
3D space simply because they are nearby objects, the gas- 
deficient enhancement in the background arises from a tight 



aggregation of galaxies segregated in the four-dimensional 
position-radial velocity phase space. Up to 15 of the galax- 
ies listed in the data set in Paper I can be identified as probable 
members of this group, since they all share similar positions on 
the sky (12 h 15 m < R.A. < 12 h 30 m and +6° < Decl. < +10°). 
In addition, 12 of these objects have LOS distances between 
27 and 30 Mpc, and 8 of them (11 out of the initial 15) have 
systemic velocities between ^600-1300 kms" 1 . Certainly, the 
lack of good resolution in the radial direction prevents us for 
claiming that we have identified a true group on a sufficiently 
safe basis. Yet, the fact that one third of its potential members 
have gas deficiencies that deviate more than 2er from normalcy 
and that two of them have Hi masses less than 10% of the ex- 
pectation values for their morphology — characteristics that are 
both typical of rich cluster interiors — reinforces the impression 
that it is not a fortuitous feature. 

3. INFALL MODEL WITH REBOUND 

To model the virgocentric velocity field, we have imple- 
mented a simple point mass model for the spherical collapse 
of a zero-pressure fluid (see, e.g., Ekholm & Teerikorpi 1994, 
and references therein). However, the solution of the problem 
is not limited to the classical deduction of the velocity field of 
objects until a singularity first develops. The calculations are 
extended beyond that time until galaxies recollapse again by 
allowing that the spherical shells of matter reaching the singu- 
larity pass through themselves and reemerge. 

In the present paper, we take the, admittedly crude, point of 
view of treating galaxies as test particles moving in a constant 
point-mass potential well, which allows us to ignore the effects 
of shell crossing on their first orbit. Shell crossing is most im- 
portant inside the previously relaxed cluster body, but galaxies 
on first infall spend little time in this region because of their 
radial trajectories and high pericentric velocities. Hence, the 
gravitational acceleration they undergo during their first cross- 
ing of the cluster core is not expected to differ substantially 
from that exerted if all the mass were concentrated at the cen- 
ter. 

Under the point mass approximation, the velocity relative to 
the center of Virgo, v, for mass-less shells of radius r in a uni- 
verse with null cosmological constant (a non-null cosmological 
constant has a negligible contribution in the local universe) is 
given by the pair of parametric equations 

77 -sin 77 _ -y/GAfyicto ^_ r sin 77(77 -sin 77) 
(I-COS77) 3 / 2 r 3 / 2 ' to (I-COS77) 2 

for the bound shell case, and by 

sinh 77-77 ^/GMyicto r smhrjismhri — 77) 

(cosh??- 1) 3 / 2 r 3 / 2 ' to (cosh?7-l) 2 

for the unbound case. In these expressions, M V ic is the effec- 
tive total mass of the VIC region (strictly speaking, the point 
mass representative of the region under consideration around 
the VIC center), to is the current age of the universe, and G 
is the gravitational constant. In our double-infall model bound 
shells reach maximum expansion for development angles 77 = it 
and 3ir, whereas full collapse is achieved when 77 = 2tt and 4tt. 

The predicted radial velocity, V, of a galaxy at a distance R 
from the Local Group (LG) and observed at an angular distance 
9 from the Virgo center can be evaluated from 
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V(R) = Vvic ± vy/ 1 -R 2 mc sm 2 0/r 2 , (3) 

where /? V ic and VVic are, respectively, the barycentric distance 
and velocity of the VIC in the LG reference frame. The lin- 
ear and angular distances of the galaxy are related through the 
cosine law 

d 2 = l+D 2 -2Dcos6 , (4) 

where d = r//?vic and D = R/Ryic are the linear distances ex- 
pressed in Virgo's distance units. The (-)-sign applies for 
galaxies with D < cos 8, while the (+)-sign is for D > cos 9. 

4. MODEL PREDICTIONS VS. OBSERVATIONS 

As shown in the previous section, our simple dynamical 
model describing the velocity field in the VIC region involves 
parameterization in terms of to, M V ic, ^vic, and Vyic — the sys- 
temic velocity can be replaced by the cosmological velocity of 
Virgo for a given infall velocity of the LG. In spite of the fact 
that only three of these parameters are independent, it is ad- 
visable to work with no more than two free parameters, as the 
current velocity-distance data on the VIC are still insufficient 
to constrain models with so much freedom. For this reason, we 
allow the Virgo mass and distance to vary freely while keeping 
Vvic fixed to 980 km s" 1 . We have, nonetheless, checked that 
the range of values of Vvic allowed by the observations does 
not lead to significantly different results. The viability of the 
solutions will be cross-checked by comparing the value of to 
predicted by the best fits to the observations with the cosmo- 
logical age of Oil Gyr derived for the £l m = 0.3, Q\ = 0.7 
CDM model in the Ho-Key Project by Freedman et al. (2001) 
and the more precise to = 13.6 ±0.2 Gyr recently inferred by 
Sievers et al. (2002) from CMB observations. 

The systemic velocity-distance diagram for the VIC region 
plotted in Figure 2a shows the basic expected features: an ini- 
tial steeply rising velocity-distance relation at the cluster front, 
a central broad region with the maximum observed velocity am- 
plitudes, and a final ascending part of the relation, expected to 
approach the local Hubble law asymptotically. Ultimately, all 
the available good-quality observations should be considered to 
define the constraints of dynamical models like ours. Unfortu- 
nately, as illustrated in Figure 2b, even after we have extended 
the predictions to the first orbit following rebound, the theo- 
retical pattern cannot explain the motions of various systems 
within the post-rebound expansion regime. This fact and the 
difficulties inherent to the modeling of the motions in the inner- 
most cluster region — where multiple rebounds are expected to 
occur — prompted us to restrict the acceptable range of mod- 
els by putting all the weight of the fits on the envelope to the 
streaming velocities within the asymptotic regime (cf. Fig. 2b). 
This is indeed the fitting procedure adopted by first-infall mod- 
els and has the advantages of: (1) relying on the subset of ob- 
servations that offer the highest constraining power; (2) being 
unaffected by shell crossing; and (3) being barely sensitive to 
the angular distance from the VIC center. 

The three solid curves in Figure 2a demonstrate the locus of 
the optimum model for lines-of-sight corresponding to angular 
separations of 4°, 6°, and 8°. We find it compelling that, in 
spite of the fact that the minimization relies on a fraction of the 
data, our best solution explains a large number of the galaxy 
motions around the VIC region. The good accordance between 



the observations and the prediction of the double-infall model is 
reinforced by the remarkable symmetry of the motions with re- 
spect to the local Hubble flow defined by the imaginary straight 
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FIG. 2. — (a) Systemic velocity vs. LOS distance for Virgo spirals. The 
curves demonstrate the predicted velocities for 9 = 4° (thick solid), 9 = 6° (thin 
solid), and 9 = 8° (dotted). The dashed portion at large Virgocentric distances 
is for unbound shells. Galaxies with uncertain distances (open symbols) have 
been excluded from the fit. Parameter DEF measures HI deficiency in units 
of the mean standard deviation for field objects (= 0.24). (b) Barycentric ve- 
locity vs. barycentric distance for spirals with reliable LOS distances. Open 
and filled symbols are, respectively, for galaxies closer and more distant than 
the barycenter. Squares identify the galaxies used to define the two asymptotic 
branches, triangles are for members of the Hi-deficient background group, and 
crosses for galaxies with DEF > 3cr and LOS distance < 15 Mpc. Error bars 
in distance are shown for the last two subsets of galaxies and for one of the 
galaxies within the asymptotic regime. They are representative of the uncer- 
tainty on the distances of all the available good-quality observations. Errors in 
the velocities are negligible in comparison. The curves are the same solutions 
depicted in (a). 

line passing through the position of the LG and that of the pre- 
dicted VIC barycenter (cf. Fig. 2b). Moreover, the error con- 
tours in the (M V ic,^vic) plane drawn in Figure 3 show that our 
predictions are in reasonable agreement with previous estimates 
of the values of these parameters. Thus, while mass calcula- 
tions of the central cluster region based on the X-ray emission 
or the virial theorem amount to about several times 10 14 M©, 
our best estimate M V ic = 2.8 x 10 15 M Q approaches the val- 
ues <2x 10 15 M Q inferred from modelings of the velocity 
field of the Local Supercluster (Tully & Shaya 1984; Fouque 
et al. 2001). The acceptable barycentric distances are also well 
within the very poorly constrained range (~ 16-24 Mpc) of VIC 
distances reported in the literature, although our most likely 
value Ryic = 21.0 Mpc advocates a large-distance scale which, 
for a cosmological VIC velocity of 1200 kms -1 , brings the lo- 
cal value of the Hubble constant to about 60 km/s/Mpc. Fur- 
thermore, our best solution leads to to = 13.5 Gyr, in excellent 
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agreement with the expansion ages advocated in Freedman et 
al. (2001) and Sievers et al. (2002). Finally, note from eq. (3) 
that the constraint provided by the motion of the LG with re- 
spect to the VIC (i.e., a null systemic velocity at our position) 
is automatically satisfied by the models. 
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FIG. 3. — 1,2, and 3cr confidence limits on Myic and Ryic obtained from 
our dynamical model of the VIC region. The cross indicates the best estimates 
of these parameters. 

5. IMPLICATIONS ON THE ORIGIN OF THE GAS DEFICIENCY 

One of the most striking predictions of our model, illus- 
trated in Figure 2b, is that various Hi-deficient galaxies in the 
frontside and backside of Virgo can be accommodated in the 
same way on an orbit following rebound and on a first infall 
trajectory. We conclude that these galaxies might not be re- 
cent arrivals but have already plunged into the VIC center in the 
past. Especially interesting is the case of the Hi-deficient back- 
ground group observed at R <~ 28 Mpc and V <~ 1000 kms" 1 , 
whose near turnaround position in the Hubble diagram and the 
assumed standard harmonic oscillation movement of the galax- 
ies around the cluster center, together with a cosmic age pre- 
sumably close to 13.5 Gyr, indicate that it might have passed 
through the Virgo core about 4.5 Gyr ago. The fit to the ob- 
served velocity field also indicates that a significant number of 
galaxies in the VIC region might be currently falling towards 
the cluster from the frontside — perhaps not for the first time 
either — , so it does not substantiate former claims of a paucity 
of objects of this sort (Fouque et al. 2001). 

We do not claim, however, that this is the only possible so- 
lution satisfying the observations. For instance, an alternative 



model based on the assumption that the gas-deficient group in 
the background is also on first infall gives similarly reasonable 
best values of the assumed parameters: M V ic = 1-8 x 10 15 M Q 
and Ryic = 20.2 Mpc, leading to to = 14.0 Gyr. Nor are we en- 
dorsing the argument that all Hi-deficient galaxies in the VIC 
region have passed through the cluster core. It is clear from 
Figures 2a,b that the positions of some gas-poor galaxies, espe- 
cially in the front cluster side, are best explained according to 
our model if they are on first infall. In any event, our dynam- 
ical modeling of the VIC region has shown that characteristics 
such as a large Virgocentric distance or a near turnaround po- 
sition are not by themselves conclusive indications of a recent 
arrival. The substantial Hi deficiency of the background sub- 
clump found in Paper I may well have originated in an earlier 
passage of this entity through the Virgo core. 

Even if our suggestion that the Hi-deficient group on the 
backside of the VIC might not be a recent arrival is finally 
proven well-founded, it is still necessary to find an explana- 
tion for the apparently long time (^4.5 Gyr) that these galax- 
ies have maintained a substantial dearth of gas without notice- 
able consequences on their morphologies: 9 of its 15 probable 
members are late-type spirals, whereas the 5 with the largest 
gaseous deficiencies have types Sb or later. (Nor do the magni- 
tudes of these latter galaxies show any evidence of a substantial 
dimming.) Yet, the details and chronology of the evolution of 
galactic properties triggered by the sweeping of the atomic hy- 
drogen, as well as its repercussions on the star formation rate, 
are still poorly understood. For instance, estimates for gas con- 
sumption time scales in the absence of gas replenishment from 
a sample of 36 spiral galaxies of various morphologies by Lar- 
son, Tinsley, & Caldwell (1980) produced values ranging from 
0.9 to 15 Gyr, with a median of 3.9 Gyr. These authors also 
discussed the color evolution of disk galaxies whose star for- 
mation has been truncated at various past times and concluded 
that substantial reddening in the latest spiral types requires that 
most of their star formation ceased about 5 Gyr ago. On the 
other hand, recent observational studies of the galaxy popula- 
tions in clusters at different redshifts (e.g., Fasano et al. 2000) 
indicate that while the quenching of star formation induced by 
the removal of the gas appears to be rapid (<~ 1 Gyr), the mor- 
phological evolution of disk galaxies takes several billion years. 
Perhaps evolution is still slower for objects, such as the mem- 
bers of the gas-deficient cloud detected in the background of 
Virgo, that spent most of their time out of the aggressive cluster 
environment. 
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